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^ ■ ABSTRACT 

o; 

Q , We present for the first time the detection of mid- J CO hne emission in the 

^ \ outskirts of an evolved massive star, which indicates the presence of warm and/or 

k> \ high density molecular gas. Aiming to learn about the interplay between evolved 

H \ massive stars and their surroundings, we have carried out CO J = 2 ^ 1, 

J = 3 ^ 2, and ^'^CO J = 2 ^ 1 line observations in a 4' x 4' field around 
the ring nebula G79.29+0.46, which is illuminated by a strong candidate to LBV 
star. The whole field shows extended predominant emission in both CO and ^^CO 
J = 2 ^ 1 lines, which probably comes from the large cloud which contains the 
star forming region DR15. When this large-scale emission is removed, minor 
scales features become evident, particularly in the CO J = 3 ^- 2 line, strikingly 
coincident with the ring nebula. The high critical density of CO J = 3 ^ 2 
(some 10^ cm~^) gives additional support for the association with the massive 
star, since high density molecular gas has more chances to survive in such a harsh 
environment. This structure may have been produced by the cumulative effect 
of a strong steady wind in the LBV stage or earlier, which has compressed the 
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surviving parent molecular cloud. In addition, immersed within this CO feature, 
we have also discovered a higher density clump (at least several ~ 10^ cm~^), 
unresolved by the telescope and probably having a higher kinetic temperature. 
Towards the clump, the presence of a jump of 14 — 16 km s~^ in the gas velocity 
may indicate the existence of a shockfront. This clump may have been created by 
at least one mass eruption, 10^ — 10"^ yr ago. Thus, this work shows that not all 
the molecular gas is destroyed during massive star evolution, and consequently 
we are dealing with a new laboratory where one can learn about the mass-loss 
phenomena associated to the brief LBV stage. 

Subject headings: circumstellar matter — ISM: general — ISM: individual 
(G79.29+0.46) — ISM: molecules — ISM: structure — stars: mass loss 



1. Introduction 



Massive stars are the most influent sources in the dynamical and chemical evolution 
of the galactic interstellar medium (ISM). Although relatively few in number, massive stars 
are the most powerful, the fastest in the evolutionary paths and the richest from a chemical 
viewpoint. Mass-loss phenomena (strong and dense winds, violent mass ejections) and a 
very intense UV field provide a significant input to their surroundings, which heats, ionizes, 
dissociates, blows out, shocks and enriches the ISM as a whole. 



According to current evolutionary models (see, for example, iMaeder &: Meynetl Il994j : 
Langer et al.lll994j ). the main sequence of massive stars lasts several 10^ yr, while the Wolf- 
Rayet stage takes some 10^ yr. All intermediate stages last typically from 10^ to several 
10^ yr. Stellar winds, size, effective temperature and UV flux vary during each phase by 
orders of magnitude, while the stellar luminosity remains almost constant. When these 
changes occur, a number of different structures in the circumstellar /interstellar gas sh ould 
be produced. Numerical simulations (JGarcia-Segura et al.lll996l : Ivan Marie et al.ll2004j ) pre- 
dict the formation of several shockfronts in the ISM at scales from tenths to tens of pc. 
The possibility of mass eruptions at the int ermediate and uiistable stages may also add a 
significant contribution to the nebular mass (jSmith &: Owockill2006l ). 



So far, the study of the ISM surrounding massive stars constitutes an excellent method 
to know about the evolutive history of the exciting source, because it is expected that 
different evolutive stages would put their "fingerprints" in the circumstellar material. Both 
neutral and ionized gas have been widely observed in the last decades, and today we can 



recognize many interstellar bubbles associated to the Main Sequence (jArnal et al.l Il999 
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Cichowolski et al.l l2003h and Wolf-Rayet nebulae flMiller fc Chul Il993l : Moore et al.ll2000h . 
among other galactic structures. In recent years, the hot dust surrounding massive stars was 
the subject of several near- and mid-IR observational campaigns, and we kn ow about the 



presence of disks, torus an d shells around stars at intermediate evolutive stages (JBecklin et al. 
199ilO'Haraet al.ll2003h . 



The search for molecular gas in these environments is still limited, possibly due to the 
a-priori hypothesis by which the parent molecular cloud around an evolved massive star 
may have been destroyed soon after star formation. Even s o, the detected low-excitation 
CO (JMarston et al.lll999l : iNota et al.ll2002l : iRizzo et al.ll2003al ) is giving indications that the 
molecular component would represent an important contribution to the total mass of gas. 
Unfortunately, the low-J lines observed to date trace densities as low as few 10^ cm~^, and 
the link with the stars or the circumstellar material can only be based on the morphology of 
the CO structures, and not on any other physical parameter. 

A special and motivating case is NGC2359. To date, this is the only Wolf-Rayet nebula 
where H2 has been detected at 2;um (ISt-Louis et al.l 119981 ). Later, t he global distributi on of 
CO has been almost si multaneously map ped using the Kitt Peak (JRizzo et al.ll2001al ) and 
SEST radiotelescopes (ICappa et al.ll200ll ). showing an important amount of molecular gas 
surroundin g the optical nebu l a. Foll ow-up studies have reported the detection of complex 
molecules (Rizzo et al.ll2001bl . l2003bl ) and multiple layers of rather dense gas (~ 10^ cm~^; 
Rizzo et al.l l2003d ). The origin of the molecular gas seems related to stages previous to 
Wolf-Rayet (dynamical lifetimes tdyn ~ 10^ — 10^ yr). 

This significant result encourages the search for molecular gas in the environs of evolved 
massive stars which have not yet become Wolf-Rayet. Of those, the luminous blue vari- 
able (LBV) stars are probably the best candidates, due to the "instant aneous" effect of their 
strong winds and spectacular mass ejections over the star surroundings (JHumphreys fc Davidson 
I994J : iNota fc Clampinl Il997l ) . The nested shells observed in several cases around LBV neb- 
ulae span masses from 0.1 to 20 Mq, most of wh ich come from short-liv ed eruptions, such as 
the well known 19th century outburst of 77 Car (ISmith &: Owockill2006l ). The multiple layer 
structure discovered in NGC2359, referred above, may have been formed during a previous 
LBV stage. 

A recent ce nsus provides a to tal of 12 confirmed galactic LBV stars and 23 candidates 
of such objects ( IClark et al.ll2005l ). G 79. 29-1-0. 46 appears as a strong LBV candidate in the 
sample. Noticed for the first time by IWendker et al.l (Il99ll ) in a radio map (1420 MHz) of 
the Cygnus X region, this ringlike structure looks highly symmetric, and has a di a meter 
of ~ 3'2. In a pioneer paper based on continuum and Hi observations, iHiggs et al.l (119941 ) 
have shown the thermal nature of the continuum emission, and suggested that the ringlike 
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nebula is an ionized shell of swept-up interstellar material; in addition, a possible Hi bubble 
is reported at a systemic velocity (LSR) of ~ —3 km s~^. By assuming a normal gas to 
dust ratio, the atomic mass is not enough to account for the extinction measured, and hence 
a significant am ount of molecular gas is predicted. After examining the IRAS High 
Resolution images. IWaters et al.l (Il996l ) concluded that G79.29+0.46 is a detached shell due 
to an epoch of high mass loss (~ 5 lO~^M0yr~^), followed by a more quiet period. They 
have estimated a shell mass of ~14M0, and measured a dust temperature of 60-77K. 
This picture is consistent with the LBV nature of G79.29+0.46's exciting star, where some 
different events of mass-loss may have occurred in the recent past. 

The central star of G79.2 9+0.46 has been spectroscopically stu died at the optical and 
the near-infrared wavelengths (jWaters et al.lll996l : IVoors et al.ll2000l ). Its LBV classification 
is supported by a high luminosity (> 10^ Lq) and a currently high mass-loss rate (> 10^^ M© 
yr~^). Its galactic longitude is not adequate to estimate a rehable kinematic distance; even so, 
Voors et al.l (120001) derived a value of fti 2 kp c. This is close to 1.7 kpc, the estimated distance 
to Cygnus 0B2 (JMassey fc Thompsonlll99ll ). a huge OB association which G79.29+0.46 may 
belong to. G79. 29-1-0.46 is behind the Great Cygnus Rift, projected 9' to the n orth-west of 
the s tar-forming region DR15, which causes high visual extinction (up to 16 mag; Kno dlsederl 
20001 ) . and where an important amount of diffuse CO were reported (ISchneider et al.. .20061 . 
and references therein). 



The aim of this work is to detect and to make a first morphological and dynamical 
analysis of the molecular gas associated to G79.29+0.46. Below we present the results of 
surveys in the rotational lines of CO J = 2 ^ 1 and 3 — i> 2, and ^^CO J = 2 ^ 1 
surrounding the LBV ring nebula. The obtained results encourage to make further studies 
about dynamics, excitation and the possible relationship between the molecular gas detected 
and the ring nebula. 



2. Observations 



The observations were made with the Submillimeter Telescope (SMT), located on Mount 
Graham, Arizona. The telescope is a 10.4 m diameter primary with a nutating secondary. 
Design, opt i cs and structural aspects of the telescope have been extensively described by 
Baars et al.l (119991 ). The single-polarization receiver SIS-230 was used to observe the CO 
and ^^CO J = 2 ^ 1 line emission; on the other hand, the CO J = 3 ^ 2 line was observed 
using the dual-channel SIS-345 receiver. The three lines were mapped in several ~ 4'x4' fields 
around G79. 29-1-0. 46, using on-the-fly mode, with a fixed position (empty of line emission) 
as the reference. Several on-the-fly maps were later averaged in order to get a final rms as 
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low as possible. The scanning directions and starting positions were different each time to 
reduce stripping. The Table 1 summarizes some parameters of the observations; for each 
line, the table indicates the frequency, angular resolution (HPBW), main-beam efficiency 
and final rms of the maps. 

The receivers were tuned for double-sideband operation, assuming a sideband ratio 
of 1 for the system temperature calibration. The receivers were calibrated with ambient- 
temperature and liquid nitrogen-cooled absorbing loads to determine the receiver noise tem- 
perature; calibration of the final data have uncertainties below 20 %. These measurements 
were made each time a receiver was tuned and at intervals of a few hours thereafter. Sys- 
tem temperatures wer e measured with the standard sky/ambient temperature load method 



(iKutner fc Ulichlll98ll ). Calibrations were performed every 15-20 minutes to track variations 
in atmospheric opacity. Typical single-sideband system temperatures during the observations 
were in the range 250 — 400 K for the 1.3 mm lines and 1200 — 2000 K for the 0.8 mm line, de- 
pending on the elevation and atmospheric opacity. Two acousto-optic spectrometers (AOS) 
have been used as backends. The AOS have 2048 channels, and a frequency resolution of 
480 kHz, providing a velocity resolution of 0.6 (0.4) km s~^ in the 1.3 mm (0.8 mm) bands, 
and bandwidths wide enough for galactic emission. The pointing was checked within 2 hr on 
a planet or other bright continuum source using the integrated line intensity. Corrections 
to the expected positions for the pointing model were not more than 4" in either azimuth 
or elevation. Planet observations were also used to check the main beam efficiency at the 
observed frequencies, which are also listed in Table 1. Throughout all the paper, the tem- 
peratures are on a main-beam scale (Tmb), and the velocities refer to the Local Standard of 
Rest system (Vlsr)- 



3. Results: emission maps around G79. 29-1-0.46 

Figure 1 shows a panel of the emission maps corresponding to the CO J = 2 ^ 1 line. 
To provide a first-look morphological comparison with G79.29+0.46, the circle in the c enter 



of each map roughly traces the external border of the ring nebula (JHiggs et al.l Il994j ) . At 
first sight, we note that large-scale emission dominates all the observed field, mostly at the 
2.4 km s~^ map. Despite of this, some minor-scale features can be observed, such as the 
unresolved arc at the southwest, around (—80", —110"). This arc is present in a wide range 
of velocities, and bounds part of the ring nebula. The inner part of the circle is filled by 
some low-level emission, especially in the range -3 to +5 km s~^. The emission at the first 
and last maps do not appear morphologically linked to the ring nebula. 

Figure 2 sketches the ^'^CO J = 2 —>■ 1 emission in a similar fashion to Fig. 1. With the 
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aim of reducing the map noise, a 2-D banning filtering has been applied to build the figure. 
This filtering slightly reduces the angular resolution of the map down to 45" . In general, the 
emission of this isotopomer is well correlated with the most intense CO emission in Fig. 1. 
Similarly, besides the most intense emission at the southeastern corner of the observed field, 
minor scale feature are also noticed. A systematic deviation of the contours is visible toward 
the north, north-east, and south-west of the plotted circle. Another outstanding feature of 
Fig. 2 is the lack of emission in the first, second, and last maps; as this behaviour can not 
be explained merely by the sensitivity of the map, it clearly indicates that the CO opacity 
is relatively higher at the central velocities (i.e., in the range —3 to +10 km s^^), where the 
isotopic ratio [CO]/[^'^CO] is lower, when compared with the more extreme velocities plotted 
in Fig. 2. 

The emission maps of the CO J = 3 ^ 2 line are shown in Fig. 3. Due to a higher 
velocity resolution, the twelve maps of the figure span the same velocity range as those of 
Figs. 1 and 2. With the exception of the last map (11.4 km s~^), the emission concentrates 
on the central velocities, mainly on the six central maps. Some extended emission appears 
at -0.5 and 1.2 km s~^, while in the following three maps the CO J = 3 — ;► 2 emission is 
concentrated outside the ring nebula. Especially outstanding is again the southwest feature; 
it has the highest level of the maps (above 10 K on a Tmb scale) , remains unresolved at 345 
GHz, and is widely extended in velocity. 



4. Linking the CO to G79. 29-1-0.46: distance, morphology, kinematics 

The field of G79. 29-1-0. 46 is immersed into Cyg X, a major star-formi ng region which 
contains large amounts of molecular gas having densities of some 10'^ cm~^ (jSchneider et al. 



20061 ). More specifically, G79.29+0.46 is projected close in the sky (only 9') to the smaller 



star-forming region DR15, although it is not clea r if they are physically connected. DR15 



is located at 1.7±0.3kpc fjSchneider et al.ll2006l ). The distance to G79.29+0.46, however 



is hard to obtain with similar precision due to the uncerta i nties in the stellar parameters 



dependent on its classification. In this sense, IVoors et al.l (J2000l ) has made an important 
contribution, providing a distance in the range 1 to 4kpc, with a likely value of 2kpc, close 
to the distance of DR 15. 

Thus, the CO emission presumably linked to G79.29+0.46 are expected to be heavily 
contaminated by foreground/background emission from DR15 or CygX as a whole. More- 
over, contribution from the solar neighbourhood does not have to be disregarded (mostly 
in the J = 2 — > 1 lines), due to the velocities involved, close to cero. As we have seen in 
Figs. 1, 2, and 3, most of the CO emission is dominated by large-scale emission which fills all 
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the mapped area. From a visual inspection to the low resolution maps of lOka et al.l (l200ll ). 
we can relate our large-scale CO emission -which increases towards the southeast- with the 
molecular cloud arou nd DR15. Furthermore, in the wide field MSX and SCUBA images of 



Redman et al.l (120031 ) we can clearly associate the CO emission to the higher column densi- 
ties of dust. In principle, we can therefore think about the large-scale emission detected in 
our maps as arising from the foreground/background emission associated to the molecular 
cloud where DR 15 is embedded. 

It is also expected that the CO gas associated to the ring nebula should have minor-scale 
features, and be morphologically connected to the ring nebula. Thus, in order to propose a 
link between the observed CO and G79. 29-1-0. 46, we should look at the smaller-scale features. 
In this sense, the southwestern clump is the most prominent feature, visible in all the maps. 
Other minor scale features can be inferred by the local deviations of the contours already 
referred to in Sect. 3. In order to better distinguish the local features presumably associated 
to the ring nebula, we have removed a plane emission (i.e., a polynomial of order 1 in both 
celestial coordinates) for each velocity channel, in the three observed lines. These planes 
have been fitted by a least-square method to the whole field, taking the external zone of 
the field as input, and assuming that every feature outside ~ 160" is purely associated to 
the foreground/background emission. The results are presented in Fig. 4, where integrated 
velocity maps in two different velocity ranges (-1-1.3 to -1-4.8 km s~^ at panels 4a to 4c; 
—3.0 to —1.3 km s~^ at panels 4d to 4f) are shown for each line. Both the positive and 
negative velocity ranges have been selected to better show the minor-scale features of the 
CO emission. The molecular emission is shown by contours, superposed on an ISO 25 /xm- 



image in grayscale (JTrams et al.l 119971 ). 



The distribution of the background-substracted emission is now closely following the 
ring-shaped geometry of the infrared nebula. At the positive velocity range (Figs. 4a to 4c) 
the most prominent feature is the southwest clump, clearly visible in both CO lines. The CO 
J = 3 ^ 2 line emission also has lower level contour emission bordering the external part of 
the nebula. The lack of a ^^CO counterpart of the southwestern clump may be due to the 
background removal method, because the intensity of the local feature may be comparable 
to the fluctuations of the fltted large-scale plane. Strikingly, at negative velocities (Figs. 4d 
to 4f), the emission of the three lines closely match the infrared nebula, especially in the 
western half of the ring. The peak of the emission is again at the southwest, but at an inner 
position with respect to the southwestern clump at positive velocities. 

The morphological coincidence of the small-scale features of the CO emission with the 
infrared ring nebula encourages their mutual association. At the best of our knowledge, the 
CO J = 1 ^ and J = 2 ^ 1 observations presented in AG Car is the only LBV nebula 



where a molecular counterpart has been reported (iNota et al.ll2002l ). In that paper, however, 
no maps were constructed and it is hard to determine if the morphology of the CO resembles 
or not that found in our data. Moreover, the low- J lines may contain important amounts of 
diffuse, background or foreground gas. In G79.29+0.46, the small-scale features that more 
clearly are correlated to the ring nebula are those of the J = 3 — > 2 line (Fig. 4c and 4f), 
and they are also the most intense. If we also take into account the higher critical density 
of the J = 3 — *> 2 line (some 10^ cm~'^), compared to the low- J lines, we can think about a 
scenario where high density molecular gas has more chances to survive in harsh environments, 
dominated by strong and variable stellar winds and UV fields from the massive progenitor 
star. 

The kinematics of the local gas was studied by position-velocity maps across selected 
strips. On many of them, we could not detect any sign of expansion or rotation of the CO 
gas. The dominant velocities across the field are from —2 to km s~^, with almost no 
systematic variations or gradients. This is, on the other hand, the velo city range where most 



of the low- resolution emission was reported as arising from DR15 (jOka et al.ll200ll ). An 
outstanding exception to this general trend is shown in Fig. 5. This figure was constructed 
by selecting the strip indicated in Fig. 4c, along a northeast to southwest direction (exactly 
45° from north to east). The slice was selected to pass through the southwest arc (at ~ 260" 
in Fig. 5), where a sudden broadening occurs. This broadening grows up to 14 km s~^ (from 
—2 to +12 km s~^), and is limited to an angular extension comparable to the beamsize. It is 
also noticeable that the CO velocity shifts when the strip passes from the inner to the outer 
region of the nebula, i.e., when passes through the southwest arc. 

The feature displayed in Fig. 5 is a dynamical indicator of a low-velocity shock, acting 
beyond the ionized ring nebula at a projected velocity of 14 km s~^. This assumption is 
supported by two facts: firstly, the region affected by the broadening is relatively thin, 
because it is unresolved by our beam; and secondly, there is no smooth variations of the 
velocity between the inner and the outer regions. An important, additional argument is 
provided by the fact that there is no CO emissio n arising from PR 15 or other close objects. 



above 10 km s~^; in this sense, the Fig. 4c of lOka et al.l (120011 ) is ellocuent, showing no 
features at velocities above 10 km s~^. In other words, the CO velocity outside the ring 
nebula is locally exceptional, and it may only be explained by some kind of interaction with 
the LBV-candidate star. The presence of shocked gas close to the hot star encourages its 
association and explains its own existence, because the (relatively) recent mass eruptions 
and variable winds after the main sequence may have produced some sort of features like 
the ones recognized in the molecular circumstellar gas. 

As were mentioned in the introduction, an increasing number of hydrodynamical sim- 
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ulatio ns are predicting the pres e nce of several shocks in the environs of evolved massive 
stars (JGarcia-Segura et al.l Il996l : Ivan Marie et al.ll2004l ). This shocked region is certainly 
of particular interest, because it may constitute a good laboratory to test these theoretical 
pr edictions. This is es pecially true when we take into account the multiple layers reported 
by lRizzo et al.l (l2003d ) surrounding a Wolf-Rayet nebula (NGC 2359), where a previous LBV 
stage was assigned as the responsible of such layers. The chemistry of the molecular gas in 
the southwestern clump may also be of particular interest, with a scenario characterized by 
a shocked cloud, and close to a highly-obscured, hot, and evolved star. 



The velocity range of the unshocked CO (from —3 to +1 km s~^) is -within the un- 
certainties of the observations and the expected dispersion- coincident with the Hi hole's 
systemic velocity reported by iHiggs et al.l (Il994j ). This Hi counterpart may represent the 
main sequen ce signature of the star progenitor, as in many other cases o f the Galaxy (see, 
for example, lArnal et al.l 1 19991 : iRizzo et al.ll2001al : Cichowolski et al.ll2003l ). T hese velocities 
do no t differ too much to those of the metal lines (IVoors et al.ll2000l ) and Bra ( IWaters et al. 
19961 ) from the star B 



Below we discuss the particular physical conditions of the molecular gas linked to 
G79. 29-1-0. 46, and also provide some quantitave arguments about its possible origin. 



5. Discussion 



5.1. Physical parameters: the LVG approach 

As we mentioned in the earlier sub-section, the detected molecular gas should have 
physical conditions which are different to those of the general ISM. The detected CO may 
simply be interstellar gas accumulated by the stellar wind, the results of several ejecta, or 
the relics of the parent molecular cloud, affected by the destructive action of UV radiation 
and/or pushing wind. An estimate of the molecular mass, together with the neutral atomic 
and the ionized gas masses, may complete the energetic balance related to the mechanical 
input from the exciting star. Other local parameters, such as the gas temperature or the 
density, will help to better describe the physical conditions in the studied region. 

In a wid e variety of gal actic scenarios, CO was found to be close to LTE but having a 
high opacity (JLequeuxl 120051 ). From the most common values of the CO to ^^CO J = 2 ^ 1 



^Note that, in both papers, the conversion from hehocentric to LSR reference systems were made in the 
wrong sense, and consequently the star and nebular velocities result to be ~ 34 km s^^ higher than the 
published ones. 
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line ratio (between 3 and 5 when ^'^CO is detected) we confirm that both CO fines are 
opticafiy tfiick. Furtfiermore, tfie fiigfi intensity of tfie CO J = 3 —> 2 fine (comparable or 
even fiigfier tfian tfie J = 2 ^ 1 fine), indicates tfiat we are deafing witfi densities (n(H2)) 
above 10^ cm~^, and rotational temperatures (T^ot) of 40 K or more. 

We can make a first approacfi to tfie pfiysical conditions of tfie CO by using a LVG code 
to predict tfie observed intensities and fine ratios. In order to compare tfie LVG modefiing 
witfi tfie observations, we fiave used tfie line-integrated maps afieady sfiown in Fig. 4 in tfie 
two velocity intervals. Tfie CO J = 3 —> 2 map fias been convolved to tfie J = 2 — > 1 angular 
resolution to work witfi a directly comparable Tub scale. Tfie CO to ^^CO J = 2 — >• 1 line 
ratio varies from 2 to 8, wfiile tfie CO J=3^2to2^1 line ratio (fiereafter referred to 
as Rexc) varies from 0.3 to 1.7 on a Tmb scale in botfi velocity intervals. 

After analyzing tfie emission in all tfie lines, as well as tfie line ratios, we fiave proceeded 
witfi tfie LVG modelling by cfioosing cfiaracteristic regions witfiin tfie observed field. Tfiese 
regions are indicated in Fig. 6, and tfie meanings are as follows: 

Region A: It contains tfie most intense emission from tfie soutfiwestern clump in tfie (1.3, 
4.8) km s^^ velocity interval. Tfiis region may fiave representative values for botfi tfie total 
mass and n{B.2) of tfie emission at positive velocities. 

Region B: Tfiis region corresponds to tfie maximum R^xc in botfi velocity intervals, and 
points out tfie maximum densities expected in tfie wfiole area. 

Region C: Tfiis region comprises most of tfie emission at negative velocities (excluding region 
B), coincident witfi tfie ring nebula. It may represent tfie most common pfiysical conditions 
in tfie velocity interval (-3, -1.3) km s^^, and may contain most of tfie mass. 

Region D: Tfiis is a control region, witfi no significant features associated to tfie ring nebula. 
Tfiis is built just to compare tfie results of LVG modelling witfi tfie otfier regions, and to 
cfieck tfie validity of tfie metfiod. 

Tfie four regions defined are sketcfied in Fig. 6 superposed to tfie distribution of i?cxc in 
tfie (1.3, 4.8) km s~^ velocity interval in greyscale, and to some contours of tfie CO J = 3 ^ 2 
emission in botfi velocity intervals (see caption). Note tfiat greyscale peaks at region B and 
not at region A, wfiere tfie line intensity is fiigfier for positive velocities. Qualitatively, region 
B migfit fiave a fiigfier density but lower mass tfian tfie region A, wfiere we fiave found tfie 
CO absolute peak and a fiigfi (but not tfie fiigfiest) i?exc- 

We fiave modelled tfie CO emission for different values of tfie kinetic temperature (Tk), 
namely 40, 70, and 100 K. For a given Tk, we fiave used tfie ^^CO J = 2 ^ 1 line intensity 
and Rcxc as input, and tfius predicted tfie ^^CO column density (A^(^^CO)) and n(H2). 
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According to the method, N(^^CO ) is most sensitive to the hne intensity, and does not vary 
too much with T^ (lLequeuxll2005l ). The obtained values for A^(^^CO) and n(H2) were later 
introduced to obtain the line opacities and excitation temperatures. The CO column density 
(■/V(CO)) was obtained from iV(^^CO) by assuming a CO to ^^CO relative abundance of 70 
(jWilson fc Matteuccilll992l ). Finally, the total molecular mass was derived from A^(CO), by 
assuming a CO abundance of 8 10~^, a distance of 1.7 kpc, and including the correction for 
the presence of He and other elements. The results from the LVG modelling are summarized 
in Table 2. For each velocity interval, we have indicated the region, its angular size, the 
assumed Tk, and the obtained values for n(B.2), A^(CO), and the total molecular mass. 
Finally, the last column provided a brief description of each region. 



The visual a bsorption Ay mea sured towards G79.29+0.46 varies from 12 (JTrams et al. 



19991 ) to 15 mag (IVoors et al.ll2000l ). By assuming a crude gas-to-dust ratio of 



A^(HI) + 2 A^(H2) ~ lO^Mv 



< 6 - 7.5 10^1 cm-2^ 



we derive A^(H2 

iV(CO) value of < 5 - 6 10^^ 



For a normal CO abundance (8 10^^), a general 
cm~^ is expected. The values of A^(CO) obtained in our 
control regions are one order of magnitude lower, which may indicate that most of the local 
gas is in atomic phase (both neutral and ionized). Moreover, the densities of the control 
regions are the lowest of the table, which again reinforces the proposed link of the molecular 
gas characterized by the regions A, B, and C. While the densities of region D are around 
10^ cm~^, the values for the other regions show and enhancement of one or two orders of 
magnitude. 

The total mass at positive velocities rounds 8Mq, roughly distributed between clumps 
A (5M0) and B (3Mq). At negative velocities, a total mass of 5Mq is computed, dis- 
tributed between the clump B (1.5 M©) and the northwestern arc (3.5 Mq). The missing 
mass of the method, mostly due to the reduction of the velocity intervals and to the selected 
regions, is in any case lo wer than 2M(7). These masses fall within the range of those ex- 
pected in L BV outbursts (ISmith &: Owockil 120061). and to other verified mass eruptions hke 
in 7] Carina (jSmith et al.l |2003| ) , or P Cyg (JMeaburnI 1200 ll ) . Although the lack of symmetry 
of the associated CO is puzzling, the hypothesis of mass ejections as the origin of such struc- 
tures is worth to be considered in the future. A follow-up study should be carried out to 
observe other emission lines at higher frequencies or larger telescopes, capable of resolving 
the southwestern structure. The dominant physical processes may also be inferred by getting 
information from molecular tracers of shocks and PDRs. 

The assumed Tk is surely not uniform along the whole field: probably, the emission 
at negative velocities is better described by higher Tk than the regions emitting at positive 
velocities. As expected, the column densities are almost independent from Tk, and the 
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variations are well below the observational uncertainties. The ^^CO J = 2 ^ 1 line was 
found to be optically thin in all cases. However, both CO lines are optically thick; the 
J = 3 ^ 2 line systematically has a larger opacity than the J = 2 ^ 1 line, which 
may result in larger values of i?exc; and consequently n(H2). The ^^CO J = 2 — > 1 line 
is subthermally excited in all cases. For positive velocities, the LVG modelling needs an 
excitation temperature above 40 K, indicating a T^ well above such value. In fact, fitting 
was impossible in region B for T^ = 40 K. 



5.2. On the origin of the detected molecular gas 

The high densities found in the field, particularly in the regions A and B, strengthen the 
hypothesis of several shockfronts acting upon the circumstellar medium around G79.29+0.46. 
A set of shockfronts is the natural consequence of the presence of di fferent wind regimes 



during star evolution, and the probable occurence of mass eruptions (JGarcia-Segura et al. 



19961 : Ivan Marie et al.ll2004l ). So far, it is important to associate the different features with 



specific evolutive stages of the star progenitor. 

We can deal with a rough estimate of the age of the high-density gas, by taking into 
account the morphology and kinematics presented above. For an assumed distance of 1.7 
kpc, the southwestern clump is projected at a distance of ~ 1 pc. A lower limit for the shock 
velocity would be 14 km s~^, the magnitude of the "jump" in Fig. 5. By adopting these 
values, a crude estimate of the dynamical time would be 

R 1 pc y, 

tdyn < T^— = ../^ , ::. 7 X 10^ yr. 
v^shock 14 km s 

This value of tdyn should be regarded as a loose upper limit for the age of the shock, mainly 
due to the fact that: {i) the present shock velocity may be higher than 14 km s~^, because 
this is the most conservative value, and just a radial projection of the real shock velocity; 
[ii) the shock might have started with values of several hundreds of km s~^, and then have 
been decelerated when travelling into the circumstellar gas. If this were the case, the real age 
would be one order of magnitude lower. By looking at the high regularity of the torus-like 
nebula, a reasonable value for the age of the structure may be below 10^ yr, perhaps some 
10^ yr. 

The youth of this structure allows us to associate the molecular shockfront with the 
actual LBV phase, which might have experienced several episodes of mass ejection in the 
near past. As the material is indeed more disperse than, for example, r/ Carina, we can 
speculate about two possibilities: {i) G79.29+0.46 is older than r] Carina, and hence has 
had more time to disperse its surroundings; or (m) the energy involved in G79. 29-1-0. 46 is 
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greater than in 77 Carina. The first hypothesis seems more plausible, because of the size of 
the circumstellar material associated to G79.29+0.46, and because recent, violent eruptions 
of rj Carina have been recorded in the 19th century. Even so, it is interesting to explore both 
possibilities, as well as to go deep into a comparison of both objects, in order to shed some 
light on the interplay between these type of stars and their outskirts. 



6. Conclusions 

A region of 4' x 4' in the field of the ring nebula G79.29+0.46 has been mapped in the 
CO J = 2 — > 1 and 3 — > 2, and the ^^CO J = 2 ^ 1 emission lines, with a maximum angular 
resolution of 24", equivalent to ~ 0.2 pc at an assumed distance of 1.7 kpc. We report for the 
first time the detection of a mid- J CO emission line from the surroundings of a LBV nebula. 
Provided the excitation conditions of the CO J = 3 — >• 2 line, the studied gas is presumably 
warm (more than 33 K) and dense (above some 10^ cm"^). As a first result of this work, 
it is shown that not all the molecular gas is destroyed during massive star evolution, and 
consequently we are dealing with a new laboratory where one can learn about the mass-loss 
phenomena associated to the brief LBV stage. 

The CO emission is dominated by large scale features, which uniformly increases towards 
the southeast; this background/foreground emission is probably associated to DR15, due to 
its close location and the velocity range of emission. Minor scale features, however, show a 
morphology completely different, dominated by a large arc which coincides with more than 
half of the infrared ring nebula, and bounding it at a projected distance of up to Ipc. This 
emission is particularly intense (in the J = 3 — > 2 line) towards the southwest, where a more 
conspicuous zone is embedded. Moreover, there is a shockfront of a projected velocity of 14 
km s~^ acting precisely in this southwestern clump. 

We have predicted the observed line emission by LVG modelling. A minimum T^ of 
40 K was determined in an extended area which have densities of the order of 10^ cm~^. 
Immersed within this plateau, there are two outstanding clumps (regions A and B), which 
have densities well above 10^ cm~^ and probably a higher T\^. 

The total molecular mass associated to the nebula is ~ 13 to 15 Mq, while the clump 
masses are about 5 Mq each. Most of the gas may result from the surviving parent molecular 
cloud, later compressed by steady stellar winds and radiation pressure. The origin of the 
clumps may be different, and we can not roule out the existence of recent mass eruptions. 
The higher density, the presence of a shockfront and the mass involved may be indicators of 
violent events, produced ~ 10^ yr ago, or less. This fact, together with the high symmetry 
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of the nebula, points to an LBV origin of the feature. 

This study of the molecular component directly associated to the LBV stage may help to 
understand the recent past of G79.29+0.46, and therefore the evolutive mechanisms of this 
kind of objects. So far, it is worth making more observational efforts in order to resolve the 
emission, particulary in the clump. The differences between regions A and B are puzzling, 
and may indicate the existence of more than one mass eruption in the recent past of the 
progenitor of G79.29+0.46. A follow-up study of complex molecules should also be carried 
out, because they can provide information about the physical parameters in detail. The 
chemistry of the molecular gas in this region may also be of particular interest, with a 
scenario characterized by a shocked cloud, and close to a highly-obscured, hot, and evolved 
star. An abundance study of such molecules may point to the dominant emission mechanisms 
in the zone, which can be related to the ionized and atomic neutral gas, together with the 
hot circumstellar dust. Finally, other surveys of CO emission in LBVs and other evolved 
massive stars will also help to understand how frequent the mass eruptions are in this class 
of objects, and how important they are when compared to the steady high mass loss and to 
the UV destructive radiation. 

We are grateful to the technical staff of SMT for their kind assistance during both the 
remote and on-site observations. We also wish to acknowledge the anonymous referee by 
his/her useful comments, which greatly improved the paper. Data reduction and analysis 
were done using the GILDAS package (http://www.iram.fr/IRAMFR/GILDAS). This work 
is supported by the Spanish Ministerio de Ciencia y Tecnica, through grant AYA2003-06473. 
The SMT is operated by the Arizona Radio Observatory (ARO), Steward Observatory, 
University of Arizona. 
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Table 1. Observational parameters 



line 






frequency 
(GHz) 


HPBW 

(") 


??MB 


rms^ 
(K) 


CO J = 


2- 


-^ 1 


230.5380000 


36 


0.85 


0.4 


i^CO J = 


2- 


-^ 1 


220.3986765 


38 


0.86 


0.2 


CO J = 


3- 


-^2 


345.7959899 


24 


0.60 


0.6 



^Approximate rms of the final maps, in T^ scale 
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Table 2. Physical parameters in selected regions 



Region Area T^ 'n{}i2) N{CO) Mass Description 
(ster) (K) (cm-''^) (cm^^) (Mq) 

Positive velocities: from 1.3 to 4.8 km s~^: 

A 5.5 10-s 40 2.4 10^ 1.4 10^^ 5.2 Southwestern clump. 

70 1.7 10^ 1.2 10^^ 4.5 CO J = 3 ^ 2 peak. 

100 9.1 10^ 1.2 10^^ 4.5 



B 3.9 10' 



D 7.6 10" 



40 


>4 10^ 


8.3 10^^ 


2.7 


i?cxc peak. 


70 


1.3 10^ 


9.5 10^^ 


3.1 




100 


2.0 10^ 


8.2 10^^ 


2.7 




40 


7.7 10^ 


5.1 10^^ 


2.5 


"Control" region 


70 


4.1 10^ 


5.0 10^^ 


2.5 




100 


2.5 10^ 


4.9 10^^ 


2.4 





Negative velocities: from -3.0 to -1.3 km s : 

B 3.9 10-s 40 2.7 10^ 4.7 10^^ 1.6 CO J = 3 ^ 2 and 

70 8.8 10=^ 4.5 10^^ 1.5 i?exc peaks. 

100 6.4 10^ 4.6 10^^ 1.5 

C 4.0 10-^ 40 1.7 10^ 1.3 10^^ 3.4 Most extended 

70 6.5 10^ 1.3 10^^ 3.4 emission, coincident 

100 4.8 10^ 1.3 10^^ 3.5 with ring nebula. 

D 7.6 10-s 40 6.4 10^ 2.7 10^^ 1.3 "Control" region. 

70 3.2 10^ 2.6 10^^ 1.2 

100 2.2 10^ 2.6 10^^ 1.2 
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Fig. 1. — CO J = 2 ^ 1 emission maps in the field of G79.29+0.46. Central velocities (in 
km s~^) are indicated on top of each map, and correspond to the mean of two consecutive 
channels. Starting and spacing contours are 0.6 K (in Tmb scale). Greyscale (in K) and beam 
size are indicated in the top left map. The central circle, ~ 4' in diameter, roughly sketches 
the external border of the ring nebula. Equatorial coordinates (J2000.0) are offsets from the 
exciting star of G79.29+0.46 
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Fig. 2. — The same as for Fig. 1, but for the ^'^CO J = 2 ^ 1 hne. Starting and spacing 
contours are 0.25 K 
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Fig. 3. — The same as for Figs. 1 and 2, but for the CO J 
spaced by 1.5 K, starting in 3K 



3^2 hne. Contours are 
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Fig. 4. — Background-removed maps of the three hnes observed, integrated in two different 
velocity ranges: from 1.3 to 4.8 km s~^ (panels a to c), and from —3.0 to —1.3 km s^^ 
(panels d to f). Lines are indica ted at the top left of each map. Contours are overlaid to an 
ISO 25 /zm image of the nebula (JTrams et al.lll997l ). Contours start and spacing are 1.6, 1.0, 
4.4, 1.0, 0.4, and 1.0 Kkm s~^ from a to f, respectively. The diagonal line on map c marks 
the strip selected for the position-velocity diagram shown in Fig. 5 
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Fig. 5. — Position- velocity diagram of the CO J = 3 ^ 2 line emission in the field of 
G79.29+0.46. The direction of the slice is traced by the line in Fig. 4c, oriented 45° from 
northeast to southwest. The notable broadening near 260" corresponds to the southwest arc 
just beyond the nebula (see text) 



-26- 



o 

CD 

o 

u 

<1 



200 


Till 


III 


II 


1 1 1 r 


100 


- 







1 




- 


D 




vOV 





— 






\yA\ 


100 


3 11 




o ^ 


1 1 1 1 1 1 1 1 r 


200 


1 1 1 1 


III 



200 100 -100-200 
Aa (arcsec) 

Fig. 6. — Regions selected for the LVG modelling, indicated by letters A, B, C and D. 
Greyscale corresponds to the CO J = 3^2to2^1 line ratio (-Rcxc) in the velocity 
interval (+1.3, +4.8) km s~^. Black, thin contours correspond to the CO J = 3 ^ 2 line 
emission in the same velocity interval. Grey, heavy contours trace the same but in the 
velocity interval (-3, -1.3) km s~^. All data have been convolved to 36" of angular resolution 



